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Abstract. We present a study of the strong coupling between radiation and matter, 
considering a system of two quantum dots, which are in mutual interaction and 
interacting with a single mode of light confined in a semiconductor nanocavity. We take 
into account dissipative mechanisms such as the escape of the cavity photons, decay of 
the quantum dot excitons by spontaneous emission, and independent exciton pumping. 
It is shown that the mutual interaction between the dots can be measured off-resonance, 
only if the strong coupling condition is reached. Using the quantum regression theorem, 
a reasonable definition of the dynamical coupling regimes is introduced in terms of the 
complex Rabi frequency. Finally, the emission spectrum for relevant conditions is 
presented and compared with the above definition, demonstrating that the interaction 
between the excitons does not affect the Strong Coupling. 
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1. Introduction 

Cavity Quantum Electrodynamics (CQED) has provided an appropriate framework to 
understand the interaction between light and matter in a full quantum level. One of 
the most relevant achievements of CQED is the coherent and reversible exchange of 
energy between an active medium and the cavity mode. This behavior is known as 
the Strong Coupling (SC) regime pQ. In high finesse QED cavities SC was achieved 
with Rydberg atoms several years ago [2, [3]. It was also realized in the last decade in 
semiconductor systems, where a possible physical system consists of a quantum dot (QD) 
embedded in a semiconductor micro(nano)-cavity [U EJ EJ [7] , in which the QD discrete 
level structure resembles the atomic CQED physics. The possibility of achieving SC in 
semiconductor systems, allows to consider possible applications such as single photon 
sources [El [9], coherence and entanglement control [TP] , quantum computation and 
quantum information processing [111 [Gl [13] , Bose-Einstein condensation of polaritons 
PU E5] or polariton lasing pm El [18] . 

Weak (WC) or strong (SC) coupling regimes for a QD microcavity system can be 
determined in a micro-photoluminescence (//PL) experiment, in which the emission 
spectrum of the micro(nano)-cavity is associated to the dynamical regimes of the system 
[T0| IT9] . In these experiments, two peaks are clearly identified far off resonance. One of 
them is associated to the excitonic transition whereas the second one is related to the 
cavity photons. Furthermore, near resonance the state structure becomes more complex 
and the simple exciton-photon picture ( WC) may change to the so-called polaritonic 
states. The splitting between the two peaks depends on the dissipation rates (decay and 
pumping), the coupling constants between the subsystems, and the detuning between 
the exciton and photon energies. Data obtained experimentally from the //PL spectrum 
can be studied as a function of the detuning parameter A. In the resonance condition 
(A = 0) if both peaks cross, then the system is in weak coupling, otherwise the system 
is in the SC regime [U [9]. 

A model that is able to reproduce the above mentioned experimental facts is presented 
in [19]. In this work, the //PL spectrum of a QD micro (nano)-cavity system is modeled 
considering photonic and excitonic incoherent pumpings, and decay processes. In 
addition, different coupling regime conditions in the linear (low power excitation density) 
and non-linear (high power excitation density) were introduced by Tejedor and coworkers 
in [20, El] , recently a similar model extended for N independent QDs coupled to a single 
common cavity mode (Dicke model) was presented in [22J. Besides these results, in 
[231 [24] the //PL spectrum associated to the coupling of two semiconductor QDs to a 
single nanocavity mode was presented. 

When more than one QD is considered they not only interact with the light mode, but 
also interact among themselves. The physics of this interaction could be used in quantum 
logic devices and quantum computation applications [25]. One of the mechanisms in 
which two excitons can interact is through a resonant energy transfer, the so-called 
Forster interaction [26J. It has has been characterized in the experiments performed in 
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[27J EE] • The aim of this paper is to show how the coupling regimes change due to the 
mutual interaction between the QDs, considering dissipative effects. Our main finding 
is that the mutual interaction between the QDs can be associated to a SC condition 
out of resonance. Nonetheless, we show that the WC and SC regimes do not change 
significantly, as a function of the mutual interaction strength between the QDs. This 
result opens the possibility of designing solid state all-optical quantum networks by 
deterministically growing QDs in nanocavities. 

This paper is organized as follows: in section 2, we present the Hamiltonian and 
dissipative dynamics of the system. In section 3, we introduce the quantum regression 
theorem (QRT), which allows to compute two-time correlation functions and the 
emission spectrum of the system. Once the theoretical tools are introduced, in section 4 
we go back to the Hamiltonian of the system, and obtain the energies and the polaritonic 
states (dressed states). Then, we show the contributions of each sub-system to the 
polaritonic states. Next, by including the dissipative effects in the dynamics of the 
system and using the QRT matrix, we introduce the complex (half) Rabi frequency 
which serves as a criterion for distinguishing the different coupling regimes. Finally to 
support our results we present the emission spectrum of the system. The discussion and 
conclusions of this work are presented in section 5. 

2. The System 

The system considered here, is composed of two interacting and spatially separated 
QDs coupled to a single nanocavity mode of a high Q photonic crystal. We model 
the QDs as two level systems, and use a Tavis-Cummings like model. Furthermore, we 
include a Forster type interaction between the QDs. This allow us to write the following 
Hamiltonian (in units in which H — 1): 

2 

H s = w a f a + ^ {cox + 9% + <?ia*\ } +9i2 + ai&t) 

The first term corresponds to the free field Hamiltonian. The creation a) and 
annihilation a operators, are associated with photons of energy u$ = ux — A (where A 
is the detuning of each exciton respect to the field mode). The first term on the sum 
is related to the exciton energy. The operators erj = |Xj) (Gi\ and a, = (JQ| are 
the creation and annihilation operators of the zth exciton. The exciton is modeled as a 
two level system, where the ground state corresponds to the absence of excitations and 
is denoted by \G), whereas the presence of an excitation in the QD i.e the exciton will 
be denoted as \X). The transition energy between these states in any QD is ux- The 
second term in the sum describes the dipolar interaction between the QDs and the light 
mode [29] in the rotating wave approximation (RWA). The strength of such interaction 
is given by Finally, the last term accounts for the Forster exciton-exciton interaction, 
with coupling constant g^- This interaction represents the resonant exchange of energy 
between the QDs. 
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On the other hand, the system-reservoir interaction Hamiltonian accounts for the next 
processes: (i) The direct coupling between each exciton and the dispersive photonic 
modes, this process is responsible for the spontaneous emission, (ii) The escape of the 
cavity mode photons, the so called coherent emission, (iii) A common continuous and 
incoherent pumping of each of the excitons. The reservoir Hamiltonian Hr contains the 
features of the environment which can be modeled as a set of harmonic oscillators. It is 
not necessary to write them explicitly. For a detailed discussion refer to the literature of 
open quantum systems [30J . Finally to make the problem tractable it is assumed that the 
interaction between the system and the reservoir is weak, and that the dynamics of the 
reservoir is memoryless, i.e the Born-Markov approximation. Under these conditions, 
the evolution of the reduced density operator of the system p$ can be written as: 



d A 

dt ps = l 



Ps,H s 



+ - [2ap s cJ - a^aps - ps^a) 



A 

+ \ Yl ( 2a *Ps & l - Oi&iPs ~ PsOi&i) (2) 



=1 

2 



where k and 7 are the coherent and spontaneous emission rates respectively, and P is 
related to the incoherent pumping rate of each exciton. The above parameters will be 
called dissipative parameters through the whole text. 

The dynamics of the populations and coherences of the density matrix can be obtained 
from equation (jSJ). For the results to be presented below we use the bare states basis: 
{\k\) ® ^2) ® \ n )}i which we will denote for convenience as: {\ki,k 2 ,n)}. In this basis 
\ki) represents the excited \X) or ground \G) states of the i-th exciton, whereas |n) 
indicates a Fock state with n photons. In figure Q] we show an schematic representation 
of the dynamics of the system. In figure 1(a) we portray all the possible transitions 
between the states of the system up to the first excitation manifold, associated to the 
Hamiltonian given by equation ([1]). It is seen that the Hamiltonian dynamics only admits 
horizontal transitions between states of the same excitation manifold, and there is no 
coupling between states belonging to different excitation manifolds. On the other hand, 
when we take into account the dissipative terms, see figured] (b), the non Hamiltonian 
dynamics couples states between different excitation manifolds. Notice that, if we 
consider only the terms k and 7, the stationary solution will always be the vacuum 
state \G, G, 0). Nonetheless, if we introduce the common pumping term the system can 
have as stationary solution populations of excitons or photons different from zero. This 
is a crucial fact in the determination of the SC, as can be seen from the results of [19|. 
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Figure 1. Dynamics of a system of two QDs in mutual interaction and interacting 
with a field mode, in presence of dissipative and pumping mechanisms. Considering up 
to one excitation, (a) The Hamiltonian dynamics Hs- (b) Non-Hamiltonian dynamics. 

3. First order correlation function and spectrum 

As we mentioned in section [TJ the characterization of the dynamical regimes can 
be obtained from the emission spectrum of the system. To obtain the spectral 
function we use the Wiener- Khint chine theorem. The emission spectrum of the 
cavity can be obtained from the first order correlation function of the field: S(u) oc 
lim^ooK [f™ dre iuJT GM(t,T)], where G (1) (t,r) = (tf(t)a(t + r)) is the first order 
correlation function. The dynamics of the first order correlation function can be 
obtained through the quantum regression theorem (QRT) [29], which states that once 
the evolution of a set of operators {C\(t)} of the form d T (Oj(t + r)) = £\ Aij{Oj(t + r)} 
is known, then, the two time expected values of Oj with any operator O also satisfies 
the same system of differential equations, i.e d T (Oj(t + t)0) = J2j Aij{Oj(t + r)O(t)}. 
We will use the QRT to obtain the equation of motion of G^'(t,r), and by studying 
the eigenvalues of the matrix that represents the two-time dynamics, which we call the 
QRT-matrix, we will introduce a SC criterion. To do so, we consider the evolution 
given by equation ([2]) including states of up to one excitation. This set corresponds to 
the {\G, 67,0), \G,G, 1), \X, G, 0) , \G,X, 0)} states, which are depicted in figured]. We 
select the set of operators: {a\ crj, <r|}, and calculate the dynamics of this operators as 
(O) = tr [Ops]- Finally, we find the QRT-matrix for the two time expected values. 

4. Results 

To understand the nature of the transitions and the emission peaks of the system 
we first study it without considering the influence of external reservoirs. We obtain 
the polaritonic or dressed states and energies of the Hamiltonian Hs, by consider 
9i — 92 = 9, (This ideal situation corresponds to the assumption that the excitons 
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have equal coupling strengths to the light). The energies are: 

\ = uj + A - g 12 , A± = (ujo + 9U 7T ^ ± K ) ; (3) 



where 1Z = y2g 2 + ((A + g 12 )/2) 2 , this result has an important meaning. If the initial 
state has one matter excitation then the squared modulus of the transition probability 
to a state with a single photon is: 

|(G,G,l|e-^ st |X,G,0)| 2 = \ (G,G,l\e- iHst \G,X,0)\ 2 (4) 

q 2 

1 (l-cos[27tt]), 



2TZ 2 

therefore the squared modulus of the transition probability oscillates with frequency 27?., 
this can be understood as the Rabi frequency of this system. In the case in which there 



is no mutual interaction between the excitons, 1Z reduces to TZ = y2g 2 + (^) , which 
except for a factor of two, due to the fact that we are considering two excitons, is the 
usual half Rabi frequency p] . Also notice that in addition to the bare detuning factor 
A, the term g± 2 enters as an extra detuning factor. In figure [2](a) we plot the energies 
as a function of the detuning A. The bare mode energies are shown in black lines, the 
cavity mode energy remains constant, while the QDs exciton energies change with the 
detuning, and intersect with the cavity mode energy at resonance. The energies Ao,± 
which correspond to equation (|3]), are shown in colors for two different situations. First 
when there is no interaction between the QDs i.e g 12 = (blue dotted lines), second, 
when the interaction parameter is turned on, in this case we set gu = 0.5g (red dashed 
lines). In both cases the lines associated to A± never cross each other. Nonetheless, 
when the QDs interaction is turned on, its effect is to slightly modify all the energies, 
moving upwards both A±, and moving downwards Ao, which causes a crossing with A_. 
Yet, another prominent aspect of the interacting case is that the minimum approach 
distance between the energies A± occur off resonance (green dashed arrow), unlike the 
non interacting case (green dotted arrow). The separation in energy can be calculated as: 
AE = A + — A_ = 271, which is a minimum for g\ 2 = —A . Therefore, the off-resonance 
minimum separation energy associated to the upper and lower polariton modes, can 
be seen as a witness of the interaction between the pair of QD excitons. Despite the 
demanding experimental conditions in the control of the interaction between the QDs, 
the value of this interaction constant can be at least in principle determined in a simple 
experimental way. 

On the other hand, the form of the eigenstates of H$, is given by: 

| Ac) =-j=(\G,X,0)-\X,G,0)), 

|A±) = -i- (~ {A + g 12 ± 211} \G, G, 1) + \X, G, 0) + \G, X, 0)) , (5) 



where N± is a normalization factor. It is seen that | Ao) is a maximally entangled 
exciton-exciton Bell state that is also completely separable from the light state. This 
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Figure 2. (a) Left. Eigenenergies of H$ as a function of the detuning A. The black 
continuous lines represent the photon mode energy (fixed at loq = for simplicity) and 
the exciton energy uj/g = Ao| gi2= o- In blue dotted lines the energies A± are plotted for 
the case g\2 = 0, as we said before they are above and below from the energy value Ao- 
In red dashed lines we plotted the energies when the interaction parameter is turned 
on gi2 — 0.5g. Green arrows indicate the minimum distance approach between the A± 
energies, dashed arrow for the case <?i2 = 0.5g and dotted arrow for the non interacting 
case gi2 — 0. Right. Eigenvectors associated to its eigenvalues A+ upper panel (b), 
A_ bottom panel (c). The red lines indicate the contributions of both of the QDs, and 
the black line shows the contribution of the light mode to the states of the coupled 
system, as a function of the interaction between the QDs. 

is a result of choosing the symmetric condition gi — g 2 — g. On the other hand 
the states |A±) have more complex entanglement properties, for instance, in the case 
where g = g\i the state |A_) reduces to the 3 qubit W-state [31]. In figure |2] (b), (c), 
we plot the contributions of the QDs and the cavity mode to the quantum states of 
the coupled system ii.c the polaritons) for each energy A± at resonance as a function 
of the relative interaction strength gu/g- The case A was not plotted because the 
occupations correspond to the maximally entangled exciton-exciton Bell state, and they 
do not depend on g 12 . The upper and bottom panels are related to the eigenenergies A± 
respectively. The black line corresponds to the occupation of the light mode, whereas 
the red line corresponds to the QDs occupations, which is the same for both of them 
because of the symmetry conditions imposed. From the upper panel (b), we can see 
that at gu = 0, the system has photon-like and exciton-like components, at gi2 = g we 
find a complete mixture of all three states, for greater values of gu/g the system goes 
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to an exciton-like state. In the bottom panel (c) we find the same condition as before 
for g u = 0, in this case the system goes to have an strong photon-like component, and 
lesser exciton-like component as g\ 2 j g increases. Now, we will obtain the emission 
spectrum for the system, this requires to take into account the dissipative effects. 
To do so, we derive the dynamical equations of the first order correlation functions 
v(t +r,t) = {(a f (t + r)a(t)), (a\(t + r)a(t)), {a\{t + r)a(t))}, by using the QRT. We 
obtain a linear system of the form : -r-v(i + r, t) = Av(i + r, t) with: 

( -2P - § + iu ig x ig% ^ 

A= ig x -\{2P + 1 )+iu x ig 12 , (6) 

\ m W12 -\i 2P + 7) +iux J 

where the parameters P, k, 7, . . . , are those already introduced in equations ([HE]). The 
dynamical evolution of v(t + r,t) is given by: v(t + r, t) = e^~ T v(t, t). Even for the 
symmetric case the complete solution is rather cumbersome and is not presented here. 
To obtain the criterion for the coupling regimes, we begin studying the eigenvalues ao,± 
of the matrix A, which are related to the positions f2o,± an d widths A 0j ± of the spectrum 
(iofo,± = fio,± + *Ao,±), we find: 

a = i (w + A - g 12 ) - ~ - (7) 

. ( g l2 + A \ 7P + 7 + « 
a± = 1 I h w ± it I 1 , 

where: 



fl =,/V-(r + i(^)), with r = ^±p2. (8) 

The number i? has a similarity with the i?a6z frequency, found at the beginning 
of this section, we called it half complex Rabi frequency, and it includes the dissipative 
effects. Notice that if R has a non-zero real part, the positions of the emission peaks 
£l± = dt(ia±) are different. On the other hand, if R is a pure imaginary number, the 
contributions to ct± will merely affect the width of peaks of the spectrum. Based on the 
last statement, we introduce a SC criterion in presence of dissipation as follows: 

The SC or WC regimes for a system of two excitons interacting with a confined 
light mode, is defined as the regime for which the complex Rabi frequency is a purely real 
(SC) or purely imaginary (WC) number, under the off-resonance condition A = — g 12 . 

Consequently, the system reaches the SC regime if T / g < y/2, this condition is 
represented by the region below the green dashed line in figure 3. However, to determine 
the coupling regime, usually the experimental spectra is analyzed near resonance. In 
this case the behaviour of the real or imaginary part (Q(R)) of the complex 

Rabi frequency R, screens the dynamical regime of the system. This can be seen in 
figure [3j in which we present a contour plot of the ratio \3t,(R)/Q(R)\ at resonance, as 
a function of the two free parameters: Y/g and gu/g- The predominance of $s(R) is 
represented in dark colors — > 0, a WC-like character), on the contrary 




Figure 3. Contour plot for the ratio |3?(i?)/3(i2)| between the real and imaginary 
parts of the complex Rabi frequency R at resonance, as a function of the decay rates 
r and the coupling constant 312 between the QDs. The SC regime is obtained for 
parameter values such that \$t(R)/Q(R)\ —> 00. The WC regime is reached when 
|5ft(i?)/3(.R)| —> 0. Green dashed line represent the upper limit for the SC regime. 
Orange, and blue dashed lines are used to represent some parameters for which we 
plot the emission spectrum in figure |4l 

the bright colors are associated to the predominance of 3?(i2) — > 00, 
a SC-like character). Note that, even below the green dashed line, $s(R) prevails and 
therefore in these conditions the characteristic two peaks feature of the SC regime is lost. 
Finally, to check that the criterion we have defined agrees among the different coupling 
regimes, in figure H] we plot the normalized emission spectrum for three different cases. 
In plots (a), (b) we choose variations in T/g, the dissipative parameters, whereas in (c) 
we take into account variations in gvi/g, represented along the vertical axis. For the 
case A = <? 12 = the emission spectra is presented in (a). The SC features arise for 
small values of the dissipative parameters and the splitting between the peaks is clearly 
observed. The symmetry in the peaks distribution is caused by the selected parameters, 
where both excitons are coupled equally to the cavity. On the other hand, as the 
dissipation parameter T/g increases, the peaks become broader and closer. Finally, for 
high values of T/g, the SC behavior is lost, then the system is in the WC regime. 

It is clear from equation (jHJ), that the effect of A and gu is to increase the splitting 
between the peaks, and therefore the SC condition can be reached by controlling these 
two parameters. However in figure H] (b), we plot the emission spectra at resonance and 
gu = 0.5(7 (represented by the vertical orange line plotted over figure EJ). As before, 
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(co-co )/g (co-o) )/g (co-co )/g 

Figure 4. Emission spectrum at resonance as a function of to — in arbitrary 
units, (a) Emission spectrum for both: the resonant condition and without interaction 
between the dots gi2 = 0, and the re-normalized complex Rabi frequency, (b) Spectrum 
for different dissipative values, keeping fixed g\2 — 0.5<?, which is related to change T 
over the orange line in figure [3] (c) Spectrum for different values of .912/5, in this case 
the dissipative parameters were keep constant T = 0.5g, and corresponds to the blue 
line in figure [3l 

we found the SC regime for small values of the dissipative parameters, and eventually 
the WC is reached for sufficiently high values of dissipation. Notice in figure EH (a), (b) 
that for the dissipative parameters T = 1, the two peaks feature is lost, but the system 
is in the SC regime (from the experimental point of view, this kind of broad spectra is 
deconvoluted into two components). Thus, again we can see that at resonance the effect 
of the interaction is to screen the SC regime of the system. 

Finally in figure H] (c) we keep constant the dissipative parameters (represented by 
the blue line over figure |3]). When g 12 = 0, the spectrum is completely symmetric. As 
gi2 increases, the splitting between the peaks increases as well, and one of the peaks 
decreases in intensity. This change in the intensity, can be explained by inspecting the 
behavior of the polaritonic states |A±) as a function of g\2- From figure |2] (b) we recall 
that |A_) acquires a strong photonic-like character as g^ increases. On the other hand 
|A + ) becomes more exciton-like; therefore the state |A_) decouples progressively from 
the transitions induced by the matter-matter interaction. 



Strong coupling of two interacting excitons 
5. Discussion and Conclusions 



11 



From the above results, it is interesting to notice that the off-resonance minimum separa- 
tion between the polaritonic energies, can be used as a witness of the mutual interaction 
between the QDs, i. e the Forster interaction, and in fact this coupling interaction con- 
stant can be measured directly. On the other hand, the most relevant result of this 
work is that we have established a criterion for the different dynamical regimes for a 
system of two interacting excitons symmetrically coupled to a photonic mode in a semi- 
conductor microcavity. This criterion is based on the real and imaginary parts of the 
complex Rabi frequency R. In accordance to experimental results of SC for this systems, 
it has been observed that the SC regime is reached for small values of T relative to g, 
the light matter coupling constant. This system will be in the SC regime for values of 
r/<? < y/2 (see the green dashed line in [3]) . On the contrary for values oiY/g > y/2 
we get WC independently of the interaction between the excitons. The above results 
could indicate that the experimental conditions to achieve Bose-Einstein condensation 
of polaritons (SC) or single photon sources (WC) are not strongly influenced by the ex- 
citon interaction. Finally from the expression of the complex Rabi frequency equation 
(jHJ), it is shown that the effect of the number of particles in the multiexcitonic system 
has two opposite effects. First, the mutual excitonic interaction is not favorable to the 
SC because it takes R to the limit of a purely imaginary number. In second place, the 
collective effects intensify the SC as the term proportional to g increases, in fact this 
term scales as a typical Tavis-Cummings factor y/Ng [32] where N is the number of 
excitons assuming the dots are identical, as we expected for a lineal regime. 
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